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ABSTRACT 

Recently much work in studying Gamma-Ray Burst (GRB) has been devoted to revealing the nature of 
outburst mechanism and studies of GRB afterglows. These issues have also been closely followed by the 
quest for identifying GRB progenitors. Several types of progenitors have been proposed for GRBs: the 
most promising objects seem to be collapsars, compact object binaries, mergers of compact objects with 
helium cores of evolved stars in common envelope episodes, and also the recently discussed connection 
of GRBs with supernovae. In this paper we consider the binary star progenitors of GRBs: white dwarf 
neutron star binaries (WD-NS), white dwarf black hole binaries(WD-BH), helium core neutron star 
mergers (He-NS), helium core black hole mergers (He-BH), double neutron stars (NS-NS) and neutron 
star black hole binaries (NS-BH). 

Using population synthesis methods we calculate merger rates of these binary progenitors and we 
compare them to the observed BATSE GRB rate. For the binaries considered, we also calculate the 
distribution of merger sites around host galaxies and compare them to the observed locations of GRB 
afterglows with respect to their hosts. We find that the rates of binary GRB progenitors in our standard 
model are lower than the observed GRB rates if GRBs are highly collimated. However, the uncertainty in 
the population synthesis results is too large to make this a firm conclusion. Although some observational 
signatures seem to point to collapsars as progenitors of long GRBs, we find that mergers of WD-NS, 
He-NS, He-BH, and NS-NS systems also trace the star formation regions of their host galaxies, as it is 
observed for long GRBs. 

We also speculate about possible progenitors of short-duration GRBs. For these, the most likely 
candidates are still mergers of compact objects. We find that the locations NS-NS and NS-BH mergers 
with respect to their hosts are significantly different. This may allow to distinguish between these two 
progenitor models, once current and near future missions, such as HETE-II or SWIFT, measure the 
locations of short GRBs. 

Subject headings: gamma ray bursts: progenitors - binaries: close — stars: evolution, formation, black 
holes, neutron, white dwarfs 


1. INTRODUCTION 

The last decade brought a great breakthrough in 
gamma-ray burst studies. The BATSE detectors on GRO 
have shown that GRBs are distributed isotropically on the 
sky, and that their brightness distribution is not consis¬ 
tent with a uniform source distribution in Euclidean space 
(Paciesas et al. 1999). Observations of GRB afterglows 
in X-ray, optical and radio wavelength domains (Costa et 
al. 1997; Groot et al. 1997b) led to identification of GRB 
host galaxies (Groot et al. 1997a) and measurements of 
their redshifts. This has solved the long standing prob¬ 
lem of their distance scale. While we learned that GRBs 
come from cosmological distances, there are still two major 
difficulties in understanding this phenomenon. First, we 
do not fully understand the physics of the outburst. Al¬ 
though several models have been proposed, they all have 
yet to meet some severe constraints imposed by observa¬ 
tions (i.e. releasing energies of 10 51 -10 54 ergs in timescales 
as short as lCU 2 s in the case of some GRBs). Second we 
do not know what are the astronomical objects leading to 
gamma-ray bursts, i.e. what are their progenitors. 

In recent years the black hole accretion disk model for 


GRBs has been given much attention (Fryer, Woosley & 
Hartmann 1999a; Meszaros 2000; Brown et al. 2000). Pro¬ 
genitors leading to this model include collapsars (Woosley 
1993; Paczynski 1998, MacFadyen & Woosley 1999) and 
binary mergers: helium star black hole (Fryer & Woosley 

1998) , double neutron stars (Ruffert et al. 1997; Meszaros 
& Rees 1997), black hole neutron star (Lee & Kluzniak 
1995; Kluzniak & Lee 1998) and black hole white dwarf 
systems (Fryer et al. 1999b). Also recently the connec¬ 
tion between supernovae and gamma-ray bursts received 
much attention (Paczynski 1999; Woosley 2000; Chevalier 
2000), however, there is still no clear evidence that these 
two phenomena are intrinsically correlated (Graziani et al. 

1999) . 

A good method of discerning among the binary progen¬ 
itors is to compare theoretical predictions of their merger 
site distributions around host galaxies with location of 
observed GRBs within host galaxies. Binary population 
synthesis can be of great help in addressing this question. 
One can calculate the properties of a given binary popula¬ 
tion and then place it in a galactic gravitational potential 
to trace each binary until its components merge due to 
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gravitational wave energy losses. This method has to deal 
however with a number of uncertainties, that are inherent 
in the binary population synthesis. Moreover there are 
uncertainties in what type and mass of a galaxy to use. 

The binary population synthesis method has already 
been applied to the study of compact object binaries in 
the context of GRB progenitors. However most studies 
have been concentrated only on double neutron stars and 
black hole neutron star systems. Lipunov et al. (1995) 
have used their “scenario machine” to model the popula¬ 
tion of double neutron star and black hole neutron star 
binaries in a galaxy. They calculated the expected log N- 
log S GRB distribution assuming that they are standard 
candles and compared it with the BATSE observations. 
Portegies-Zwart and Yungelson (1998) have considered the 
origin and properties of double neutron star systems and 
black hole neutron star binaries. They considered a few 
binary population synthesis models, with varying kick ve¬ 
locities, initial binary separations, initial mass ratios dis¬ 
tribution and also considered cases with and without hy¬ 
per accretion in the common envelope stage. They found 
the rates of mergers to be consistent with the GRB rate 
provided that GRBs are collimated to about ten degrees, 
and mentioned that double neutron stars may travel Mpc 
distances out of a Milky Way like galaxy before merg¬ 
ing. Bloom, Sigurdsson & Pols (1999) considered double 
neutron stars as possible GRB progenitors, and calculated 
distributions of mergers of these binaries around galaxies 
with different masses, varying the average kick velocities 
in the code. They found that a significant fraction of dou¬ 
ble neutron stars merge outside their host galaxies. Bulik, 
Belczynski & Zbijewski (1999) considered the mergers of 
binaries containing neutron stars, and Belczynski, Bulik 
& Zbijewski (2000) investigated differences between the 
populations of black hole neutron star binaries and the 
double neutron stars. Belczynski et al. (2000) found that 
black hole neutron star binaries merge closer to the hosts 
than the double neutron stars. Fryer et al. (1999a) consid¬ 
ered other types of binary progenitors of GRBs within the 
framework of the black hole accretion disk model of the 
GRB central engine. These were white dwarf black hole 
mergers, helium star black hole mergers, and collapsars in 
addition to the double neutron star systems and black hole 
neutron star binaries. They performed a thorough param¬ 
eter study, and repeated the calculations with a number of 
modifications of their standard evolutionary model. They 
calculated the distribution of merger sites in the poten¬ 
tials of galaxies with the masses of Mmw (Milky Way), 
0.25 Mmw, and 0.01 Mmw, however they do not vary the 
galactic size with mass. Bloom, Kulkarni & Djorgovski 
(2001) presented a very detailed study of the observational 
offsets between observed afterglows and GRB hosts galax¬ 
ies. They compare these observations with the theoreti¬ 
cal distributions calculated with the code of Bloom et al. 
(1999) and conclude that the so called delayed merging 
remnants i.e. double neutron star systems and black hole 
neutron star binaries are unlikely to be GRB progenitors, 
and argue in favor of the prompt bursters like collapsars 
and black hole helium star mergers. 

In this work we extend our previous studies (Belczynski 
& Bulik 1999; Bulik et al. 1999; Belczynski et al. 2000) to 
include four more proposed binary progenitors: compact 


object (black hole or neutron star) white dwarf binaries 
and Helium star mergers with black holes or neutron stars. 
We use a much improved and well tested binary popula¬ 
tion synthesis code and for consistency we also present 
the updated results for the two previously studied types 
of proposed progenitors: double neutron star and black 
hole neutron star systems. We calculate the properties of 
the ensemble of each type of the proposed GRB progeni¬ 
tors and find their distributions around different types of 
host galaxies. We compare the observed GRB distribution 
around host galaxies with the models. In order to verify 
the robustness of the results we perform a detailed param¬ 
eter study and discuss the population synthesis models 
which are responsible for the largest differences. 

An additional way of telling which group of the proposed 
binaries might be responsible for GRBs is to predict their 
rates and compare them to the observed rate of GRBs. 
Population synthesis is a powerful tool for predicting rates 
of binary populations although it suffers from many uncer¬ 
tainties as some parameters of single and binary evolution 
are poorly known. Moreover, population synthesis works 
well in predicting the relative numbers of events, while cal¬ 
culation of absolute rates requires additional assumptions. 
However, such attempts have been made by a number of 
authors mentioned above. Using the population synthesis 
method we calculate merger rates of white dwarf neutron 
star, white dwarf black hole, double neutron star, neutron 
star black hole systems, and formation rates of Helium 
star black hole and neutron star mergers. We compare the 
BATSE detection rate of GRBs with the cosmic rates of 
the binary progenitors predicted in our calculations. 

In § 2 we describe the population synthesis code Star- 
Track used to calculate properties of binary GRB progeni¬ 
tors, in § 3 we present the results, and finally § 4 is devoted 
to discussion with conclusions. 

2. POPULATION SYNTHESIS MODEL 
2.1. Stellar Evolution 

We use the StarTrack population synthesis code (Bel¬ 
czynski, Kalogera & Bulik 2001b). Here, we summarize 
only the basic assumptions and ideas of the code. 

The evolution of single stars is based on the analytic for¬ 
mulae derived by Hurley, Pols & Tout (2000). With these 
formulae we are able to calculate the evolution of stars for 
Zero Age Main Sequence (ZAMS) masses: 0.5 — 100 M© 
and for metallicities: Z = 0.0001 — 0.03. We follow 
the stellar evolution from ZAMS through different evo¬ 
lutionary phases depending on the initial (ZAMS) stel¬ 
lar mass: Main Sequence, Hertzsprung Gap, Red Giant 
Branch, Core Helium Burning, Asymptotic Giant Branch, 
and for stars with their hydrogen-rich layers stripped off: 
Helium Main Sequence, Helium Giant Branch. We end 
the evolutionary calculations at the formation of a stellar 
remnant: a white dwarf, a neutron star or a black hole. 
There are two modifications to the original Hurley et al. 
(2000) formulae concerning the treatment of (i) final rem¬ 
nant masses, and (ii) Helium-star evolution (see Belczynski 
et al. 2001b; Belczynski & Kalogera 2001). 

The StarTrack code employs Monte Carlo techniques to 
model the evolution of single and binary stars. In this 
work we use StarTrack to evolve large ensemble of stars, 
and calculate statistical properties of the binary GRB pro- 
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genitors. 

A binary system is described by four initial parame¬ 
ters: the mass M\ of the primary (the component which is 
initially more massive), the mass ratio q between the sec¬ 
ondary and the primary, the semi-major axis of the orbit 
A, and the orbital eccentricity e. Each of these initial pa¬ 
rameters is drawn from a distribution and we assume that 
these distributions are independent. More specifically, the 
mass of the primary is drawn from the Scalo initial mass 
function (Scalo 1986), 

$(Mi)«Mf 2 ' 7 (1) 

and within the mass range M\ = 5 — 100 M & . The distri¬ 
bution of the mass ratios is taken to be 

*(</) = !, 0 < <7 < 1, (2) 

following Bethe & Brown (1998). The initial binary sepa¬ 
rations are assumed as in Abt (1993) 

r 0*) oc ^ , (3) 

and finally, the initial distribution of the binary eccentric¬ 
ity is taken following Duquennoy & Mayor (1991) 

77 (e) = 2e , 0 < e < 1 . (4) 

As we are interested only in the systems capable of 
producing binary GRB progenitors, containing at least 
one NS or BH, we evolve only massive binaries, with 
primaries more massive than 5 M 0 . During the evolu¬ 
tion of every system we take into account the effects of 
wind mass-loss, asymmetric SN explosions, binary interac¬ 
tions (conservative/non-conservative mass transfers, com¬ 
mon envelope phases) on the binary orbit and the binary 
components. We also include the effects of accretion onto 
compact objects in common envelope (CE) phases (Brown 
1995; Bethe and Brown 1998; Belczynski et al. 2001b) and 
rejuvenation of binary components during mass transfer 
episodes. Once a binary consists of two stellar remnants 
(NS, BH, WD), we calculate its merger lifetime, the time 
until the components merge due to gravitational radiation 
and associated orbital decay. 

The StarTrack code may be used in several tens of 
modes, allowing for the change of main evolutionary pa¬ 
rameters and initial distributions. In the following, to¬ 
gether with the given above initial distributions, we define 
standard evolutionary model, with the set of parameters, 
thought to represent our best understanding of stellar sin¬ 
gle and binary evolution. (1) Kick velocities. Compact 
objects receive natal kicks, when they form in supernova 
explosions. Neutron star kicks are drawn from a weighted 
sum of two Maxwellian distributions with a = 175kms _1 
(80%) and cr = 700kms _1 (20%) (similar to the one of 
Cordes & Chernoff 1998). For black holes formed via par¬ 
tial fall back we use smaller kicks, but drawn from the same 
distribution as for NS. The kick scales with the amount 
of material ejected in SN explosion or inversely with the 
amount of falling back material (i.e., bigger the fall back, 
smaller the kick). For BH formed in direct collapse of 
massive stars, we do not apply any kicks, as no super¬ 
nova explosion accompanies the formation of such objects. 


(2) Maximum NS mass. We adopt a conservative value of 
A/ nl ax .ns = 3 M 0 (e.g., Kalogerafe Baym 1996). The mass 
of a compact object is estimated based on the mass and 
evolutionary status of its immediate progenitor, and not 
on any a priori assumptions. Once the mass of a compact 
object is calculated, its type (either NS or BH) is set by 
the value of M maXi ns- Thus the choice of M max ,NS does 
not affect the overall population of compact objects. How¬ 
ever, it affects the rates and various properties of binary 
GRB progenitors, as some of their groups contain either 
NS or BH (see §3). (3) Common envelope efficiency. We 
assume «ce x A = 1.0, where a is the efficiency with which 
orbital energy is used to unbind the stellar envelope (e.g., 
Webbink 1984), and A is the measure of the central con¬ 
centration of the giant (e.g., Dewi & Tauris 2000); (4) 
Non-conservative mass transfer. In cases of dynamically 
stable mass transfer between non-degenerate stars we al¬ 
low for mass and angular momentum loss from the binary 
(see Podsiadlowski, Joss, & Hsu 1992), assuming that the 
fraction / a of the mass lost from the donor is accreted to 
the companion, and the rest (1 — / a ) is lost from the sys¬ 
tem with specific angular momentum equal to 2ixjA 2 /P. 
We adopt / a = 0.5 (e.g., Meurs & van den Heuvel 1989) 
and j = 1 (e.g., Podsiadlowski et al. 1992); (5) Star for¬ 
mation history. We assume that star formation has been 
continuous in the disk of a given galaxy. To assess the 
properties of the current population of GRB progenitors, 
we start the evolution of a single star or a binary sys¬ 
tem tbirth ago, and follow it to the present time. The 
birth time fbirth is drawn randomly from the range 0-10 
Gyr, which corresponds to continuous star formation rate 
within the disk of our Galaxy (Gilmore 2001). (6) Ini¬ 
tial Binarity. We assume a binary fraction of /bi = 0.5, 
which means that for any 150 stars we evolve, we have 50 
binary systems and 50 single stars. (7) Metallicity. We 
assume solar metallicity Z = 0.02. (8) Stellar Winds. The 
single-star models we use (Hurley et al. 2000) include the 
effects of mass loss due to stellar winds. Mass loss rates 
are adopted from the literature for different evolutionary 
phases. For H-rich massive stars on MS (Nieuwenhuijzen 
& de Jager (1990); using Z dependence of Kudritzki et al. 
1989), for RG branch stars (Kudritzki & Refiners 1978), on 
AGB (Vassiliadis & Wood 1993) and for Luminous Blue 
Variables (Hurley et al. 2000). For He-rich stars W-R 
mass loss is included using the rates derived by Hamann, 
Koesterke & Wessolowski (1995) and modified by Hurley 
et al. (2000). 


2.2. Dynamical Evolution of Stars in Model Galaxies 

The population synthesis code allows us to calculate the 
age for each system at the time when both stellar remnants 
have formed, and the subsequent merger time of a given 
system based on the remnant masses and their orbit. We 
also calculate the systemic velocity gain due to asymmetric 
SN explosions and/or associated mass loss. We use these 
informations to propagate binary GRB progenitor systems 
in different galactic potentials, and to compute the distri¬ 
bution of their mergers sites around different mass and size 
hosts. 

The potential of a spiral galaxy can be described as a 
sum of three components: bulge, disk, and halo. A good 
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representation of the galactic disk and bulge potential was 
presented by Miyamoto & Nagai (1975): 


*(R,z) 


_ GMi _ 

y/R 2 + (a,. + 


( 5 ) 


where the index i refers to either bulge or disk, cq and bi 
are the parameters, M is the mass, R = y/ x 2 + y 2 and 
the x — y coordinates span the galactic plane. The dark 
matter halo potential is spherically symmetric 
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where r c is the core radius. The halo potential corresponds 
to a mass distribution p = p c /[ 1 + (r/r c ) 2 ], and we intro¬ 
duce a cutoff radius r cu t beyond which the halo density 
falls to zero, in order to make the halo mass finite and the 
halo gravitational potential is <f>(r) oc r -1 when r > r cu t- 

We consider galaxies with four masses, expressed in 
the units of Milky Way mass (Mmw = 1-5 x 1O 11 M 0 ): 
1.0, 0.1, 0.01, 0.001 x Mmw- For the Milky Way mass 
galaxy the bulge potential (i = 1) is described by: ai = 
Okpc, bi = 0.277kpc, Mi = 1.12 x 10 10 M 0 ; the disk 
potential (i = 2): a 2 = 4.2 kpc, b 2 = 0.198kpc, M 2 = 
8.78 x 1O 1O M 0 ; the halo potential: r c = 6.0kpc, and 
Mh = 5.0 x 10 10 Mq and the r cut = lOOkpc (Paczyn- 
ski 1990; Blaes & Rajagopal 1991). To obtain the poten¬ 
tial of a galaxy with the mass cxMmw we rescale all the 
masses by the factor of a and the distances a,i,bi,r c , R, z 
by a 1//3 . Such scaling keeps the galaxy density constant, 
and we made sure that our model galaxies have flat rota¬ 
tion curves. 

We adopt the following distribution of stars within the 
disk of a given galaxy (Paczynski 1990): P(R, z)dRdz = 
P(R)dRp{z)dz. The radial distribution is exponential 

P(R)dR(xRe~ R/R ^dR (7) 


and extends up to i? max - The vertical distribution is also 
exponential 

p(z)dz oc e~ z/z “ xp dz. (8) 

For a Milky Way type galaxy we have R eX p = 
4.5kpc, i? max = 20kpc, z exp = 75pc, and these parameters 
are assumed to scale with the galaxy mass as a 1 / 3 . 

Each binary moves initially with the local rotational ve¬ 
locity in its galaxy, and has no vertical component of veloc¬ 
ity. After each supernova explosion the kick imparted on 
the binary is added and the binary trajectory is calculated 
until the merger occurs. 


3. RESULTS 

3.1. Binary GRB Progenitor Types 

Fryer et al. (1999b), suggested the possibility that white 
dwarf neutron star (WD-NS) merger may lead to forma¬ 
tion of a black hole accretion disk system, followed by a 
GRB. Since the GRB outburst mechanism is not well un¬ 
derstood, the results of hydrodynamic calculations of stel¬ 
lar mergers should be treated with some caution. However, 
we will consider the group of WD-NS systems as poten¬ 
tial GRB progenitors, for the sake of completeness of the 


study. From the entire group of coalescing WD-NS bina¬ 
ries we chose these, which have the best chance to produce 
observable GRBs, i.e., systems in which WD are more mas¬ 
sive than M rn j n \y fj = 0.9M 0 to make sure that the mass 
transfer onto the NS is unstable, and the total mass of the 
system satisfies Mwd + Mns > Mmax,NS + 0.3 M 0 , since 
we require that the NS has to accrete enough material to 
collapse to a BH and that the disk formed in the merger 
must have the mass of at least ~ 0.3 M 0 to produce a 
GRB (Fryer et al. 1999b, 2001). 

Fryer et al. (1999b) also suggested that mergers of white 
dwarf black hole (WD-BH) binaries may give a rise to 
GRBs. Following the Fryer et al. (1999b) hydrodynami- 
cal calculations, we require that WD mass is larger than 
O.9M 0 to classify a coalescing WD-BH system as a poten¬ 
tial GRB progenitor. Only for these high WD masses, the 
accretion of WD onto BH is dynamically unstable and the 
rapidly disrupted (in several binary rotations) WD forms 
a thick disk around BH, which may give a rise to GRB. 

In Figure [l] we present the distribution of WD masses in 
coalescing WD-NS and WD-BH systems. For WD-BH bi¬ 
naries, WD masses distribution rises sharply at ~ O.3M 0 
and then falls approximately exponentially to flatten out 
for masses higher that ~ 0.7M 0 . Therefore, changing 
the M min yvD to slightly higher/smaller values, will de¬ 
crease/increase the number of WD-BH GRB progenitors 
roughly proportionally to the WD mass limit change. For 
WD-NS systems, the WD masses distribution is also rather 
flat close to and over M min ^wD = 0.9M 0 . However, as 
the value of limiting mass of WD is highly uncertain, we 
will present two models with the decreased and increased 

-M min t WD ■ 

Fryer & Woosley (1998) proposed yet another type of bi¬ 
nary GRB progenitors, i.e. binaries merging in CE events, 
with one component being an evolved (giant) star and the 
other already a compact object, either NS or BH. In this 
scenario, the binary does not have enough orbital energy 
to eject the common envelope, so the compact object spi¬ 
raling in finally merges with the helium core of the giant. 
The compact object disrupts tidally the helium star, ac¬ 
creting part of its material, and becoming a BH, if it was 
not one already. The remainder of the giant’s helium core 
forms a thick accretion disk around the BH, a configura¬ 
tion, which is believed to give a rise to a GRB. 

In our models we distinguish systems that contain ei¬ 
ther a NS or a BH at the onset of the CE phase, leading 
to the final merger. We will denote systems containing NS 
as helium star/neutron star mergers (He-NS), and con¬ 
taining BH as helium star/black hole mergers (He-BH). 
Following the detailed studies of He-NS and He-BH merg¬ 
ers (Bottcher & Fryer 2000; Zhang & Fryer 2000), we 
choose only these systems in which helium core mass ex¬ 
ceeds = 6 Mq as GRB progenitors. 

The distributions of helium core masses for He-NS and 
He-BH at the onset of the final CE phase are shown in Fig¬ 
ure |j. Both distributions rise sharply at ~ 1M 0 , then for 
He-NS mergers the distribution falls down rapidly above 
5 Mq, while for He-BH mergers the distribution decrease 
starts at a higher mass, 7M 0 , and is more gradual. Thus 
the number of GRB progenitor He-NS mergers depends 
strongly on the value of M m i nt H e , so we will present mod¬ 
els with different value of this limiting helium core mass. 
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Finally, the most intensively studied binary progeni¬ 
tors of GRBs are double neutron stars (NS-NS) and neu¬ 
tron star black hole systems (NS-BH). In defining the 
boundary between NS-NS systems and NS-BH systems 
we assume that the maximal mass of a neutron star is 
Af ma x,NS = 3Mg, however we will also present results for 
two smaller limiting masses of 2.0,1.5 M©. The mass dis¬ 
tribution of compact objects in these types of binaries, 
starts in our code with the maximum at the ~ 1.2 Mq, fol¬ 
lowed by a rapid decline (which reflects the shape of the 
assumed initial mass function) and then at around 3 Mq 
flattens out and stays roughly constant up to the highest 
BH masses of ~ 14 Mq. The maximum BH mass is set by 
the effect of wind mass loss on massive stars. This distri¬ 
bution is presented and discussed in detail by Belczynski 
et al. (2001b). 

Belczynski & Kalogera (2001) and Belczynski, Bulik & 
Kalogera (2001a) identified new subpopulations of NS-NS 
binaries. The new subpopulations dominate the group of 
coalescing NS-NS systems, and moreover, they were found 
to exhibit quite different properties than the systems stud¬ 
ied to date. Given the importance of these subpopulations 
to our conclusions, in the following subsection we briefly 
summarize the results of Belczynski & Kalogera (2001) 
and Belczynski et al. (2001a). 

3.2. Double neutron star binaries 

Double neutron stars are formed in various ways, includ¬ 
ing more than 14 different evolutionary channels, identified 
in Belczynski et al. (2001b). We find that the entire pop¬ 
ulation of coalescing NS-NS systems, may be divided into 
three subgroups. 

Group I consists of non-recycled NS-NS systems (con¬ 
taining two non-recycled pulsars), which finish their evo¬ 
lution in a double CE of two helium giants. Two bare CO 
cores emerge after envelope ejection, and they form two NS 
in two consecutive SN type Ic explosions. Provided that 
the system is not disrupted by SN kicks and mass loss, 
the two NS form a tight binary, with the unique charac¬ 
teristic that none of the NS had a chance to be recycled. 
For more details see Belczynski & Kalogera (2001). Group 

II includes all the systems that finished their evolution 
through a CE phase, with a helium giant donor and a NS 
companion. During the CE phase a NS accretes material 
from the envelope of the giant, becoming most probably 
a recycled pulsar. The carbon-oxygen core of the Helium 
giant forms another NS soon after the CE phase ends. The 
system has a good chance to survive even if the newly born 
NS receives a high kick because it is very tightly bound af¬ 
ter the CE episode. For more details see Belczynski et al. 
(2001a). Group III consists of all the other NS-NS systems 
formed, through the classical channels (e.g. Bhattacharya 
& van den Heuvel 1991). 

In our standard model, group II strongly dominates the 
population of coalescing NS-NS systems (81%) over group 

III (11%) and I (8%). This is due to the fact that we al¬ 
low for helium star radial evolution, and usually just prior 
to the formation of a tight (coalescing) NS-NS system we 
encounter one extra CE episode, as compared to the classi¬ 
cal channels. This has major consequences for the merger 
time distribution of the NS-NS population, and in turn 
for the distribution of NS-NS merger sites around their 


host galaxies. Merger times of classical systems are com¬ 
parable with the Hubble time, and that gives them ample 
time to escape from their host galaxies. As it has been 
shown in previous studies (e.g., Bulik et al. 1999; Bloom 
et al. 1999) which did not include Helium star detailed 
radial evolution, a significant fraction of the NS-NS popu¬ 
lation tended to merge outside host galaxies, exactly like 
group III - the classical systems. In contrast, the bina¬ 
ries of Groups I and II, due to the extra CE episode, are 
tighter, and their merger times are much shorter: of or¬ 
der of ~ 1 Myr. Thus even if they acquire high systemic 
velocities due to the asymmetric SN explosions, they will 
merge within the host galaxies, near the places they were 
born. Group I and II dominate the population, and thus 
the overall NS-NS distribution of merger sites will follow 
the distribution of primordial binaries or star formation 
regions in the host galaxies. 

The formation of NS-NS systems of group I and II de¬ 
pends on the assumption that evolved low mass helium 
donors can initiate and survive CE phase. This assump¬ 
tion have yet to be proven by detailed hydrodynamical 
calculations. 

Since, the properties of these new subpopulations were 
already discussed separately (Belczynski et al. 2001a) we 
present here only the results for the overall population of 
NS-NS binaries. 

3.3. Characteristic Binary Timescales 

We call the time a given system needs to evolve from 
ZAMS to form two stellar remnants, the evolutionary time 
and denote it by t eV oi- We call the time required for these 
stellar remnants to merge due to gravitational radiation 
the merger time and denote it by t merg . The total lifetime 
of a given system is the sum of the two: t\a e = t e voi + tmerg ■ 
In Figure 0 we show the distributions of both evolutionary 
and merger times for all GRB binary candidates, while in 
Table [I] we list the distributions medians and spans defined 
as the time range containing 90% of the systems around 
the median. 

In general, the evolutionary delays are of the order of a 
few to several tens Myrs, and their distributions are rather 
narrow for different types of systems. Since the rate of evo¬ 
lution of a given star depends primarily on its mass, t eV oi 
is set mainly by the mass of a given binary components. 
The evolution proceeds slower for less massive stars, and 
f eV oi is determined in general by the mass of the secondary 
(unless the mass ratio is reversed due to the mass trans¬ 
fer). This is why t evo i for WD-NS and WD-BH systems is 
the longest and almost equal (~ 26 Myr), as WD are the 
lightest components of binary GRB candidates. The NS- 
NS and NS-BH binaries are formed in shorter times, with 
the median of t eV oi distributions of ~ 19 and ~ 8 Myr, re¬ 
spectively. Evolutionary times for He-NS and He-BH t evo \ 
are very short (~ 9 Myr), as they finish their evolution 
even before formation of a second remnant. 

Merger times are quite different for various binary GRB 
candidates. For He-NS and He-BH mergers we do not list 
tmerg, as these events take place even before two stellar 
remnants are formed, due to the components merger in 
CE spiral in. The shortest merger times are found for NS- 
NS binaries and for WD-NS systems, with the medians of 
~ 0.7 and ~ 6.8 Myr, respectively. Much longer merger 
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times are characteristic of WD-BH systems (~ 97 Myr), 
with the longest t merg found for NS-BH binaries (~ 535 
Myr). 

3.4. Event Rates 

The method of population synthesis requires use of quite 
a number of parameters, and initial distributions of vari¬ 
ables which may affect the final results. In order to assess 
their influence on the final results we have repeated the 
calculations with varying evolutionary parameters. The 
models and the differences with the standard model are 
listed in Table ||. The coalescence rates of different types 
of GRB progenitors within each model are shown in Ta¬ 
ble |j. They have been calibrated to the Type II supernova 
empirical rates and normalized to our Galaxy (Capellaro, 
Evans & Turatto 1999). The standard model (A) results 
are based on a simulation of 3 x 10 7 binaries, while the 
remaining models are the simulations of at least 10 6 bina¬ 
ries. The statistical accuracy of most rates is better than 
a few percent, however in some cases where the rates are 
smaller than 1 Myr -1 the accuracy is of the order of a few 
tens percent, yet improving them would require a huge 
computational effort. 

Models Bl-13 represent the results of evolution with dif¬ 
ferent kick velocities imparted on the compact objects. In 
model B1 we assume symmetric SN explosions, whereas in 
models B2-12 we draw the kick velocity 14 from a single 
Maxwellian: 

g(y k ) cx 14 2 exp [—(Vfc/cr) 2 ] , (9) 

varying a values in the range 10 — 600 kms -1 . In model 
B13 we use a kick distribution of the form suggested by 
Paczynski (1990): 

/Wfc)oc [l + (V k /a) 2 }-\ (10) 

which allows for a significant fraction of low-magnitude 
kicks. We use a = 600 kms -1 , which gives a reasonable 
fit to the population of single pulsars in the solar vicinity 
(Hartman 1997). 

In Figure [| we show the dependence of coalescence rates 
on the assumed kick velocity distribution. Due to the fact 
that high kicks tend to disrupt binaries, the number of 
systems formed with at least one compact object, falls off 
quickly with the kick velocity. This has already been no¬ 
ticed for NS-NS and NS-BH binaries (e.g., Lipunov, Post- 
nov & Prokhorov 1997; Belczynski & Bulik 1999), and 
also for WD-BH systems and He-BH mergers (Fryer et 
al. 1999a). As seen from Figure [|, the coalescence rates 
of GRB progenitors fall off approximately exponentially 
with the kick velocity. Note however, that the slope is 
smaller for He-NS and He-BH mergers than for other types 
of progenitors. This is due to (i) only one kick the system 
receives, (ii) the relatively high total mass of the binary 
(recall that M rn .; n He = 6.0M 0 ), so the kick imparted to 
NS or BH does not have a big impact on such systems. 

For the majority of models the rates for WD-NS merg¬ 
ers stay close to several coalescence events per Myr per 
Galaxy. However, the rate changes significantly (0.03 — 
114.7 Myr -1 ) for a few extreme models. Besides the strong 
dependence of the rates on kick velocity discussed above, 


a large number of WD-NS is produced in models D1 and 
D2, in which the neutron star maximal mass is smaller 
than in the other models. This increase is due to our re¬ 
quirement that only systems with total mass higher than 
the maximal neutron star mass increased by 0.3 Mg are 
classified as GRB progenitors. 

In fact, the coalescence rate of all WD-NS systems (ir¬ 
respective of WD or total system mass) is as high as 204.5 
Myr -1 and 128.9 Myr -1 for our standard model and model 
D2, respectively. For the very low CE efficiency of model 
El, the rate drops down almost to zero, since in this model 
many systems potentially able to form a WD-NS GRB pro¬ 
genitor evolve through CE phase. Once the CE efficiency 
drops, a system needs to use more orbital energy to expel 
the envelope, and it becomes tighter. At very small ef¬ 
ficiencies, there is not enough orbital energy for envelope 
ejection, and the two stars merge, thus decreasing the final 
number of WD-NS systems. 

Coalescence rates of WD-BH vary much less than these 
of WD-NS progenitors, and they stay close to few coales¬ 
cence events per Myr per Galaxy for most models. The 
largest change (0.2 — 22.1 Myr -1 ) appears for models D2 
and El similarly to the case of WD-NS. The lowest rate 
of model El is explained in the same way as for WD-NS 
progenitors. The highest rate of model D2, reflects the 
fact of the lowest NS/BH mass limit, and many systems 
classified in other models as WD-NS, here are counted as 
WD-BH. 

The coalescence rates of He-NS change almost by 3 or¬ 
ders of magnitude (0.1 — 73.6 Myr -1 ), although for most 
models, including the standard one, they remain close to 
several events per Myr per Galaxy. The number of he¬ 
lium star mergers strongly (2 orders of magnitude) de¬ 
pends on the required minimum mass of the Helium core 
- see models Pl-2 and also Figure [| This is explained by 
our adopted IMF, which gives more low mass stars, there¬ 
fore low-mass helium cores are much more abundant (note 
the high rate of model PI with M m i nt H e = 4.0 Mq) than 
the massive ones (low rate of model P2 with = 

8.0Mq). The lowest rates are found for models D1 and 
D2, with the maximal NS mass of 2.0 and 1.5 Mq. With 
decreased M moXi jvs, as compared to our standard model, 
we choose only the lightest possible primaries, which will 
evolve to form NS. On the other hand we require that the 
secondary must form 6.0 Mq helium core, so it needs to 
have been massive already at the start. In these models 
only the binaries with relatively comparable mass compo¬ 
nents (q ~ 1) may evolve to form He-NS mergers. Thus, 
decreasing M max> jvs, narrows down the range of q in which 
He-NS mergers may be formed, which results in drop of 
their rate (recall that we adopted a flat initial mass ratio 
distribution). 

The rates for He-BH mergers are the most independent 
of model parameters varying just by an order of magnitude 
(4.5 — 91.6 Myr -1 ), which considering the extreme changes 
in model parameters and initial distributions is quite re¬ 
markable. The coalescence rate for standard model is 23.5 
events per Myr per Galaxy, and it remains approximately 
at this level for the majority of models. 

The dependence of NS-NS and NS-BH merger rates on 
the model parameters is discussed in detail by Belczynski 
et al. (2001b), and in what follows we restrict the descrip- 
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tion just to a brief summary. 

Merger rates for NS-NS systems change by two orders of 
magnitude in various models (2.5 — 302.2 Myr -1 ) and the 
standard model rate is about 50 merger events per Myr per 
Galaxy. As these systems experience two SN explosions, 
and the NS receive highest possible kicks (not lowered as 
in the case of BH), their rate depends very strongly on 
the assumed kick velocity distribution. The highest rates 
are found for smallest kick models (Bl-5). Production of 
coalescing NS-NS binaries is greatly reduced by reducing 
CE efficiency (model El), for the reasons described above 
(see the discussion of WD-NS merger rates). Also alter¬ 
ing the distribution of the initial mass ratio (model M2) 
changes the rates significantly and leads to an enhanced 
production of NS-NS systems. 

Finally, the merger rate of NS-BH systems stays at a 
rather constant level (1.3 — 36.2 Myr -1 ), with most model 
rates of approximately 10 events per Myr per Galaxy. The 
rate is not so sensitive to the kick velocity as the merger 
rates of NS-NS systems, because NS-BH binaries receive 
at least one smaller kick (that imparted on the BH) and 
also NS-BH systems are more massive, so the kicks have 
smaller chance to disrupt them. The smallest merger rate 
is found for the model with the enhanced wind mass loss 
rate (model G2). Due to the high mass loss the stars do 
not form massive compact objects, and the number of BH 
formed (and systems harboring BH) is greatly reduced. 
The highest rate is achieved by the shift in NS maximal 
mass of model D2 to its lowest value adopted here, which 
enhances the rate of NS-BH and depletes the rate of NS- 
NS systems. 


3.5. Redshift distribution 

Standard Model The results of the population synthe¬ 
sis code can be combined with the cosmic star formation 
rate history to yield the rate of various types of GRB pro¬ 
genitors as a function of redshift. Star formation history 
at high redshift is not well known, however it is gener¬ 
ally agreed that the star formation rate rises steeply up to 
2 ~ 1. At higher redshifts the analysis of the Hubble Deep 
Field (Madau et al. 1996) provided lower limits on the rate, 
yet these limits decrease with increasing redshift. On the 
other hand Rowan-Robinson (1999) argues that the star 
formation does not decrease and remains roughly at the 
same level above 2=1. We consider two cases: a star for¬ 
mation function falling down steeply above 2 ~ 1 (the thin 
line in Figure |^), and a case of strong star formation con¬ 
tinuing up to 2 = 10 (the thick line in Figure ||). We adopt 
a flat cosmology model with density parameter of mat¬ 
ter S7 m = 0.3, density parameter of cosmological constant 
Qa = 0.7 and for Hubble constant Ho = 65km s -1 Mpc -1 . 

For a given type i of the GRB progenitor we can calcu¬ 
late the number of events up to the redshift 2 per unit of 
observed time: 

/ , . f z 0 dr z RAz) , „. 

rateA< z) = 47t / r? —- -dz , (11) 

Jo dz 1 + 2 

where r z is the effective distance r z = cHq 1 1 + 

2 3 ) -f Ha) - 1 / 2 d 2 , and c is the speed of light. Ri{z) is the 


rate of a given type of event at the redshift of 2 : 

/•t(z=oo) 

Ri{z)= / Rsfrit') f l pi{t{z) -t')dt' , (12) 

Jt{z) 

where t is the dynamical time, dt = —Hq( 1 + 
2) -1 (fi m (l + 2 3 ) + fiA) -1//2 e?2, pi(t) is the probability den¬ 
sity of a merger of a given type as a function of time since 
formation of the system, and /) is the mass fraction of 
the binaries in the entire stellar population (single and 
binary) of mass range [0.08-100 M 0 ] that can lead to for¬ 
mation of GRB progenitors of type i. R s f r (t) is the cosmic 
star-formation rate at a time t or a corresponding redshift 
2 . We obtain the probability density Pi(t) numerically for 
each type of a merger using the population synthesis code. 
In calculation of pi ( t ) we take into account both, the evo¬ 
lutionary time delay (from formation of the system until 
two components form stellar remnants) and the merger 
time delay (the time needed for two stellar remnants to 
merge due to gravitational wave emission). 

The redshift dependence of GRB progenitor rate is pre¬ 
sented for our standard model and for the two adopted 
star formation rate histories in Figure |^. For any given 2 
the GRB progenitor merger rates are the highest for NS- 
NS binaries, than for He-BH and He-NS mergers, which 
are closely followed by the NS-BH systems. We find the 
lowest rates for mergers of WD-NS and WD-BH binaries. 

The shape of the star formation rate determines the 
shape of the GRB progenitor rate redshift distribution. 
For the Rowan-Robinson (1999) SFR, progenitors are ex¬ 
pected even at very high redshifts (2 ~ 10), while for the 
SFR of Madau et al. (1996) we do not expect to produce 
any GRBs from binary mergers over 2 <j 4. However, 
GRBs are observed at high redshifts. The highest spectro¬ 
scopic redshift 2 = 4.500 ± 0.015 was measured for GRB 
000131 (Andersen et al. 2000), while Fruchter et al. (1999) 
estimated photometrically the redshift of GRB 980329 to 
be ~ 5 (although following Bloom et al. 2001, in the Ta¬ 
ble [| we list for this burst a more moderate estimate of 
2 3.5). Therefore, if we assume that GRBs originate 

in binary progenitors our results argue against the SFR 
drawn along the lower limits of Madau et al. (1996), while 
GRBs observed at high redshifts are in agreement with our 
results based on SFR of Rowan-Robinson (1999). 

The curves in Figure |(] can be compared with the 
BATSE gamma-ray burst detection rate corrected for 
BATSE sky exposure, which is ~ 800 events per year. 
Only the rate of NS-NS and He-BH mergers is signifi¬ 
cantly above the BATSE observed rate, if we count the 
merging events up to the highest GRBs observed redshifts 
of 2 = 4 — 5. The He-NS and NS-BH merging rates in 
our standard model are marginally consistent with the ob¬ 
served rate, and these only for the Rowan-Robinson (1999) 
SRF model. Within the standard population synthesis 
model progenitors with WD merge at considerably lower 
rates than that expected for GRB progenitors. The pre¬ 
dicted cumulative rates presented in Figure^ will decrease 
if we account for collimation and thus restricted visibility 
of gamma-ray bursts. Since GRBs are thought to be colli¬ 
mated (Harrison et al. 1999; Stanek et al. 1999; Kuulkers 
et al. 2000; Panaitescu & Kumar 2001) this puts further 
limits on the binary progenitors. If any degree of collima¬ 
tion is taken into account, we may also argue against NS- 



BH and He-NS mergers besides WD-NS and WD-BH bi¬ 
naries as the sole progenitors of GRBs. Moreover, for most 
frequent NS-NS and He-BH mergers it would be difficult 
to reproduce the observed GRB rate with any significant 
degree of GRB collimation. The total rate of all binary 
mergers is ~ 7000 per year (up to z = 5), and the colli¬ 
mation which would reduce this number to the observed 
BATSE rate would be 0 w 25° (the outflow half-opening 
angle). Of course, if any single binary merger model was 
to reproduce the observed rate, the required collimation 
would be much smaller (i.e., 0 much bigger). 

Parameter Study The redshift dependence of GRB cos¬ 
mic rate is presented for all the different evolutionary mod¬ 
els listed in Table || and for Rowan-Robinson (1999) SFR 
in Figure [fj. For each progenitor type, we see that there 
are models which fail to reproduce and at the same time 
there are always models which exceed, sometimes signifi¬ 
cantly, the observed rate. Moreover the rates for most of 
progenitors are very sensitive to the assumed evolutionary 
model. Therefore, due to the population synthesis uncer¬ 
tainties, we are not able to confirm or reject any binary 
GRB progenitors just purely on the basis of their rates. A 
similar note of caution should be added to any conclusions 
about GRB collimation based on the population synthe¬ 
sis results (e.g. Lipunov et al. 1995, Portegies-Zwart and 
Yungelson 1998). The intrinsic spread in the rates when 
considering different population synthesis models is up to 
two orders of magnitude which corresponds to a factor of 
ten in the estimates for collimation. 


3.6. Distribution around host galaxies 

Standard Model In the standard model (model A) we 
have evolved N tot = 3 x 10' initial binaries and 4577 WD- 
NS, 2369 WD-BH, 9656 He-NS, 23494 He-BH, 52599 NS- 
NS and 8105 NS-BH coalescing systems formed. Next, 
we distributed the systems in a galactic disk and assigned 
galactic velocities and propagated until the merger times 
as described in § 2.2. Besides the GRB progenitor systems, 
for each galaxy mass we also propagated a number of coa¬ 
lescing WD-WD binaries, to trace the galactic distribution 
of stars. 

We show the results of the propagation calculations in 
Figures fl 10 In each figure we show the cumulative 


distributions of the projected distances of a given type 
merger. The projected distance is the distance in the di¬ 
rection perpendicular to the line of sight and we have av¬ 
eraged over all possible orientations of the host galaxy. In 
each figure, we also show, with a dashed line, the initial 
stellar distribution within a galaxy of a given mass. Note 
the different cut off radius (i? ma x) of the initial distribu¬ 
tion at 20.0, 9.3, 4.3 and 2.0 kpc for the four galaxy masses 
defined in § 2.2. 

The case of NS-BH mergers is shown in Figure |j. The 
mergers spread out with decreasing mass of the galaxy 
and even in the case of a large galaxy a significant number 
of NS-BH mergers takes place outside of the host. For a 
massive galaxy ( Mmw ) 20% of NS-BH mergers will take 
place outside the disk of the host, and as much as 70% 
will escape hosts of small mass (0.01 — 0.001 x Mmw)- 
This is due to the kicks that lead to velocities above the 
host escape velocity and relatively long lifetimes of NS-BH 
binaries. 


Mergers of WD-BH binaries take place within massive 
hosts while a significant fraction escapes from low mass 
galaxies (see Figure |j). For galaxy masses 1 — 0.1 x Mmw 
almost all WD-BH system mergers trace their initial dis¬ 
tribution. For galaxy masses of 0.01,0.001 x Mmw , 15% 
and 35% WD-BH mergers take place outside of hosts, re¬ 
spectively. These systems receive at most one kick during 
the evolution, and the gain of the velocity is not large 
enough for these binaries to escape from the potential well 
of a massive galaxy. On the other hand, WD-BH systems 
have rather long lifetimes, and if the potential well is not 
deep enough to keep them inside the galaxy, they escape 
and merge far away from the galaxies, as in the case of 
small mass hosts. 

Lighter WD-NS systems tend to merge within host 
galaxies, with only a slight dependence on the host mass 
(see Figure 0). For massive galaxies, all their mergers 
take place close to the places they were born, while for 
smallest galaxies up to 10% merge outside but close to the 
host outer regions. Since they are lighter than the WD- 
BH systems, and on average they receive higher kicks, one 
could expect that their mergers should be spread out more 
than these of WD-BH binaries. The distribution of the 
merger sites for a given mass galaxy is in general the re¬ 
sult of two competing effects; (i) the magnitude of kicks 
the systems of a given type receive and (ii) the systems 
characteristic lifetimes. These two effects are not indepen¬ 
dent; the binary lifetimes become smaller with stronger 
kicks because then only the tight, strongly bound systems 
survive. As it turns out, for WD-NS systems the short 
lifetime effect dominates over the velocity effect, and al¬ 
though they receive higher kicks they do not have enough 
time to travel outside the host before the merger takes 
place. 

Locations of He-NS and He-BH merger sites follow 
closely the initial distribution of their birth places, inde¬ 
pendently of the host galaxy mass (see Figure [Tcj) . This is 
primarily due to their very short lifetimes but also to their 
small systemic velocity gain. Both, He-NS and He-BH sys¬ 
tems have the shortest lifetimes of all the potential GRB 
progenitors studied here (see Table [lj). Their mergers take 
place even before the secondary finishes its nuclear evolu¬ 
tion (i.e., before it forms a remnant) in the CE phase, when 
the secondary evolves off the main sequence and expands 
to giant size. For all the other progenitor types, both stars 
have to first form the stellar remnants, and then usually 
considerable time is needed for gravitational radiation to 
bring the two remnants together to a final merger. Also 
He-NS and He-BH systems are relatively heavy, so the one 
kick the system experiences, does not have a great effect 
on the systemic velocity. 

Distribution of the projected distances of NS-NS merger 
sites follows very closely the initial distribution of primor¬ 
dial binaries (see Figure |l(]). Only 2% of NS-NS stars 
merge outside a massive host, and as little as 8% escape 
and merge outside of the lightest dwarf galaxies. These 
systems receive two kicks, however due to their very short 
merger times of the order ~ 1 Myr (see Belczynski et al. 
2001b for a discussion of the merger time distribution) they 
predominantly merge within even the smallest hosts. The 
NS-NS merger site distribution is quite similar to that of 
WD-NS systems. However, for NS-NS merger distribution 
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there is a tail extending to large distances from the host 
for small mass galaxies. This tail is due, to these systems 
which were formed along classical channels. These NS-NS 
have much longer merger times (typically 1 — 10 Gyr) than 
the rest of systems formed through one of the newly recog¬ 
nized pathways. The small contribution of these systems 
to the entire population of coalescing NS-NS binaries does 
not change the overall tendency of NS-NS to merge within 
even the lowest mass hosts. 

Parameter Study To study the dependence of our results 
on the assumed evolutionary parameters and initial distri¬ 
butions, we have calculated distributions of GRB progen¬ 
itor merger sites for all models listed in Table [|. For all 
our models, and for all simulated galaxy masses, He-NS 
and He-BH mergers follow the initial distribution of ini¬ 
tial binaries, and merge within their host galaxies. For 
all other systems, the results of our calculations are pre¬ 
sented in Figure [ll], [l^, for two extreme host galaxy masses 
of 1.0 Mmw and 0.001 x Mmw- 

The distribution of merger sites of WD-NS systems is 
rather independent of the model parameters and these sys¬ 
tems merge mainly within host galaxies, irrespective of the 
host mass. Most of the models are concentrated around 
the standard model distribution. Just in a few models 
more than 10-15% of WD-NS mergers take place outside 
of the smallest hosts. The two most extreme cases were 
identified in Figure [Tl] and they correspond to models N 
and F2. Model N represents nonphysical case of stellar 
evolution (and shall be treated as such), in which no he¬ 
lium giant radial evolution is allowed. This model was 
calculated just for comparison with previous results which 
did not take in to account this effect. Model F2, represents 
evolution in which every mass transfer episode (except CE 
phase) is treated conservatively, i.e., all material lost from 
the donor is accreted by the companion ( f a = 1). The 
effect of such a treatment, as compared to our standard 
evolution where half of the material is lost form the system, 
is that post-MT systems have wider separations, since no 
material and thus no angular momentum is lost from the 
binary. Naturally, the final WD-NS binaries are wider as 
well, and have longer merger times, which allow some sys¬ 
tems to escape from host galaxies. Such a model is rather 
extreme, as we know that during MT events material is 
lost from at least some systems (e.g., Meurs & van den 
Heuvel 1989). 

Distributions for WD-BH merger sites show quite signif¬ 
icant spread, allowing the possibility that majority of these 
systems merge outside the low-mass hosts. Although for 
massive hosts, most of the models show that these systems 
merge within the host boundary, for low mass galaxies as 
much as 40%, or even more merge outside the hosts. Two 
most extreme cases are these for models designated as E3 
and LI. For both models, the systems formed after CE 
or a MT phase are wider than for our standard evolution¬ 
ary scenario. Therefore, it is natural that WD-BH bina¬ 
ries have longer merger times and have greater chances 
of escaping the hosts. As we double the CE efficiency to 
acE x A = 2) in model E3, during the CE phase binaries 
use much less of their orbital energy to expel the common 
envelope. Due to this smaller energy loss post-CE bina¬ 
ries are left with wider orbital separations than they would 
have for «ce x A = 1 of our standard model. Decreased 


to half of its value the angular momentum loss, j = 0.5 
of model LI, directly influences the separations of post- 
MT systems. And although, in this model some material 
is lost from the systems during MT, unlike for model F2 
discussed above, the angular momentum loss is much de¬ 
creased, so binaries are much wider than for our standard 
model (j = 1). 

Mergers of NS-NS predominantly take place inside host 
galaxies. For massive host, all models follow very closely 
the initial distribution binaries, and depending on the 
model 95% or more of mergers take place within massive 
hosts. For lowest mass galaxies, all but two models, give 
~ 90% or more mergers within a host boundary of 2 kpc. 
Models E3 and F2 stand out, but still even for these two, 
more than ~ 80% of NS-NS mergers happen within the 
smallest mass hosts. In the above, we haven’t taken into 
account model N, marked in Figure [l^ with a dot long 
dashed line. As mentioned before, this is an nonphysical 
model, and is shown here just for comparison with previ¬ 
ous results. In agreement with previous calculations (e.g., 
Belczynski et al. 2000) for the case of massive galaxy, in 
model N about 30-50% of NS-NS stars merge outside host, 
or further away from host than 10-20 kpc. For lowest mass 
galaxies, in model N, as many as ~ 70% of NS-NS mergers 
take place outside hosts, or a few kpc from the center of 
the hosts. Detailed discussion of this significant change of 
results for NS-NS is presented in Belczynski et al. (2001a). 

Distribution of NS-BH merger sites around host galax¬ 
ies is quite sensitive to the model parameters. Although 
for the case of propagation in the potential of a very mas¬ 
sive galaxy at least 70% of these systems merge within 
hosts, majority of these mergers takes place far away from 
low-mass hosts. The curves corresponding to models D1 
and D2 clearly differ from all the remaining distributions. 
Since these two models, have lower maximum NS mass 
(ilL ma x,NS = 2.0,1.5M 0 ), many objects classified in the 
standard model (M maXj ns = 3.0M 0 ) as NS-NS are in¬ 
cluded as NS-BH in the distributions of models D1 and 
D2. These is the reason why the NS-BH distributions in 
models D1 and D2 resemble the standard model NS-NS 
distribution. If in fact, the maximum neutron star mass 
is much lower than our assumed 3.0M 0 , most of the NS- 
BH are expected to merge within even small galaxies, with 
only the heaviest binaries escaping their hosts. 

3.7. Comparison of the merger sites with GRB 
observations 

The discovery of gamma-ray burst afterglows by the 
Beppo SAX satellite have lead to the identification of GRB 
host galaxies, and to the localization of GRB events with 
respect to these galaxies. In Table ^ we list the data on 
GRB positions around host centers. Most of these are 
taken from Bloom et al. (2001), and we have added the 
entries for three recent bursts. 

From Table ^ we see that GRBs take place not far 
from the centers of their host galaxies. For some bursts, 
GRB970508, GRB000418 and GRB010222, the offsets are 
very small, and positions of the optical afterglows are co¬ 
incident with host centers. Moreover, the host galaxies 
are typically small, irregular, with intense star formation 
(e.g., Fruchter et al. 1999; Holland 2001; Ostlin et al. 2001; 
Bloom et al. 2001). One has to note that the data pre- 
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sented in Table describes only the long GRBs since only 
for these bursts afterglows have been observed so far. 

Comparing the theoretical distributions like these cal¬ 
culated above with observations is a difficult task. Ide¬ 
ally one would like to compute the theoretical distribu¬ 
tion of angular offsets between the GRB and its nearest 
galaxy, taking into account the fact that the nearest galaxy 
may not necessarily be the host galaxy. Such a calcula¬ 
tion would require a number of assumptions, about the 
evolution of galaxies with redshift, about the rate of star 
formation in galaxies, and about the mass and size distri¬ 
bution of galaxies as a function of redshift. Each of these 
quantities is uncertain by itself. Thus a calculation like 
that in our current knowledge of the evolution of galax¬ 
ies would depend on a number of uncertain assumptions 
and thus could lead to very uncertain results. However, 
results of calculations, which take into account some of 
the listed above effects, were independently obtained and 
presented by Perna & Belczynski (2001). Here we adopt 
a more straightforward approach. We assume a cosmo¬ 
logical model with Hq = 65 km s -1 Mpc -1 , Hm = 0.3, 
and Ha = 0.7, and calculate the physical distance to the 
galaxy claimed to be the host galaxy. We then use the 
Kolmogorov-Smirnov test (e.g., Press et al. 1992) to ver¬ 
ify the hypothesis that the observed distribution of offsets 
has been drawn from the distribution of offsets for a given 
type of GRB progenitor around a galaxy of a given mass. 
In each case we repeat such calculation for a number of 
population synthesis models listed in Tabic || to assess the 
range of systematic errors introduced by the population 
synthesis. 

We present the results of these calculations in Table || 
For each progenitor type we list the Kolmogorov-Smirnov 
test probabilities in case of the four galaxy masses defined 
in §2.2. We also list the highest and the lowest probabil¬ 
ity obtained when different models (with the exception of 
nonphysical model N) of population synthesis were used. 
Table g allows to evaluate the viability of each type of the 
GRB progenitor. 

Let us assume in this discussion that we reject a given 
hypothesis if the KS test probability is below 1%. One 
thing becomes immediately clear from Table || i.e. GRB 
progenitors do not reside in large galaxies like the Milky 
Way. GRB afterglows are related to small galaxies with 
the masses around 0.01 Mmw (0.015 x lO n M 0 ). This has 
been noted by the observers claiming that the typical host 
galaxy mass will lie in the range of 0.001 — 0.1 x lO n M 0 
(e.g., Ostlin et al. 2001, Bloom et al. 2001). It has to 
be noted that for each type of progenitor the value of the 
Kolmogorov-Smirnov test probabilities is a strong function 
of the model galaxy used. 

Double neutron star mergers are an acceptable choice 
and in the case of a low mass galaxy 0.01 Mmw, the prob¬ 
ability that the observed offset distribution is the same as 
the theoretical one is very high. Thus inclusion of the ad¬ 
ditional formation channels for this type of binaries has 
a significant effect (e.g., this possibility was rejected by 
Bloom et al. 2001). The range of probabilities covered by 
several different models of population synthesis is small. 
In our models the population of NS-NS mergers is domi¬ 
nated by short lived systems. Only in the case of a very 
low mass galaxy 0.001 Mmw, do the kick velocities play a 


role. Here the lowest KS test probability corresponds to 
the model B1 with very small kick velocities Bl. Since, 
for the smallest kicks, NS-NS binaries form in wider orbits 
and with longer lifetimes (e.g., Kalogera 1996) they have 
more time to escape from their host galaxies. 

The Kolmogorov-Smirnov test results presented in Ta¬ 
ble H show that we can certainly reject NS-BH mergers as 
GRB progenitors. The highest probability is obtained for 
the case of a O.IMmw mass galaxy, but its value is still 
not acceptable. This is due to the fact that these binaries 
are rather long lived (see Table [l]), therefore NS-BH can 
escape from host galaxies and merge far away from host 
centers. KS test probabilities rise to acceptable values in 
models D1 and D2, where the maximum mass of a neu¬ 
tron star is lower. In such models a number of binaries 
classified typically as double neutron stars contribute to 
the NS-BH population. 

The values of the probabilities for the mergers involving 
white dwarfs (WD-NS and WD-BH) are large and make 
these models acceptable. In the case of WD-BH mergers 
the probability even rises to 0.94 for the case of a very small 
mass 0.001 Mmw galaxy, however this number is rather 
uncertain due to the very wide range of K-S probabilities 
obtained for different population synthesis models. These 
uncertainties are not that large in the case of 0.01 Mmw 
and larger galaxies. 

Similarly we can not reject the He-BH and He-NS merg¬ 
ers. In these cases the probabilities are not as large as in 
the case of WD mergers, however these groups also con¬ 
stitutes a viable GRB progenitor. One should note that 
the systematic errors due to different population synthesis 
models are very small in this case, and the main factor 
that influences the value of the KS test probability is the 
distribution of stars in a model galaxy. He-BH and He-NS 
mergers evolve on very short timescales and therefore take 
place in star forming regions. 

3.8. Comparison with other studies 

Merger Rates A comparison of our merger rates of NS- 
NS and NS-BH binaries with a number of other studies has 
been discussed in detail by Belczynski et al. (2000b). In 
short, our rates are in good agreement with previous the¬ 
oretical predictions. Although we have noted some signif¬ 
icant differences, we attributed them to the more approx¬ 
imate treatment of stellar and binary evolution in earlier 
studies and to our recognition of new NS-NS populations. 
Recognition which was based on the assumption that CE 
phases initiated by evolved low-mass helium stars do not 
always lead to binary component mergers. 

Merger rates for several other binary GRB candidates 
have been so far presented only by Fryer et al. (1999a). We 
are not able to directly compare the rates because Fryer 
et al. (1999a) did not define the masses of WD selected to 
enter the WD-BH GRB progenitor candidate group. We 
encountered a similar problem in the case of He-BH merg¬ 
ers, for which the masses of He cores are not given. How¬ 
ever, if we assume that their systems correspond to our 
definition of GRB binary candidates, then we note a very 
close resemblance of the He-BH the rates, and a rather 
good agreement of WD-BH rates. Fryer et al. (1999a) 
found their rates of 0.15 Myr -1 for WD-BH and 14Myr _1 
for He-BH for their standard evolutionary model. In our 
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models closely resembling the Fryer et al. (1999a) standard 
model, in particular these with smaller binary fraction 
(Kl) and higher kicks (B9-12), we predict WD-BH rates 
of 0.4 — 1.0 Myr -1 and He-BH rates of 5.8 — 12.8 Myr^ 1 . 

Host Merger Site Distributions We may compare our 
results to those of Bloom et al. (1999), Bulik et al. 
(1999), Fryer et al. (1999a), and Belczynski et al. (2000). 
These authors calculated distributions of NS-NS and NS- 
BH mergers around different mass galaxies. The main 
conclusion of these studies was that a significant fraction 
(up to 40% for massive hosts, and up to 80% for low mass 
hosts) of NS-NS and NS-BH binaries merge outside of host 
galaxies. For NS-BH binaries we find very good agreement 
with previous studies, as we find that up to 25% and 80% 
of these binaries will merge outside a massive and low mass 
host (see Figure [T^). Although our calculations show big¬ 
ger concentration of NS-BH mergers in massive hosts, it is 
explained by the fact that we have adopted decreased kicks 
for BH, and therefore systemic velocity gain is decreased 
as well. However, our conclusions for NS-NS mergers are 
very different from all previous studies, due to the newly 
recognized short lived populations of these binaries, as dis¬ 
cussed throughout this work. 

Both Fryer et al. (1999a) and Bloom et al. (2001) as¬ 
sumed that He-BH mergers will take place in the star for¬ 
mation regions of the host galaxies. With our calculations 
we may confirm that, in fact, He-BH and He-NS merger 
sites, follow exactly star formation regions in their host 
galaxies. 

Fryer et al. (1999a) argued that WD-BH merger sites are 
also concentrated within host galaxies, a conclusion that 
was later adopted by Bloom et al. (2001). Our detailed 
calculations show that, in fact, for massive galaxies these 
systems follow closely the initial primordial binary distri¬ 
bution, and merge within hosts. However, for small mass 
galaxies, a significant fraction of WD-BH binaries merge 
outside of hosts. Depending on the assumed evolutionary 
model as many as 50% of these systems may merge outside 
of small mass hosts (see Figure 0 ) as discussed in §3.6. 
Fryer et al. (1999a) argued that as these systems have very 
short merger times of ~ 100 Myr, and they will not have 
enough time to escape from hosts. We find that the merger 
times of these systems are indeed of order of ~ 100 Myr 
(see Table |l|) . Nevertheless, the actual calculations of WD- 
BH trajectories prove that this conclusion is not valid for 
small mass and size hosts (0.01 — 0.001 x Mmw)- 

4. DISCUSSION AND CONCLUSIONS 

We have presented calculation of rates and spatial dis¬ 
tributions around host galaxies of several binary merger 
events, which were proposed as possible GRB progenitors. 
We have used the StarTrack population synthesis code in 
our calculations. 

We have found that the rates are very sensitive to the 
assumed set of stellar evolutionary parameters. Using the 
rates alone we were not able to exclude any of the pro¬ 
posed binaries as GRB progenitors, since the highest rates 
obtained were always higher than observed BATSE GRB 
rate for any type of a binary. In the framework of the stan¬ 
dard population synthesis model (model A) we find that 
the total rate of all the proposed binary events is roughly 
ten times larger than the observed GRB rate. However, 


we find that the spread in the rates due to uncertainties 
in population synthesis is large, and in some cases exceeds 
a factor of ~ 100. This corresponds to the uncertainty 
in the estimate of the collimation of a factor of ten. On 
the other hand we note that our standard model (model 
A) leads to an expected collimation half-opening angle of 
O ;> 25°. The measured collimation angles are somewhat 
smaller than this value, typically a few degrees (Panaitescu 
& Kumar 2001). Estimates of the GRB rates or collima¬ 
tion based on population synthesis alone carry large sys¬ 
tematic errors. 

Distributions of binary system merger sites around 
galaxies may be compared to the locations of GRB op¬ 
tical afterglows with respect to the galaxies identified as 
their hosts. Most of GRBs take place inside or close to 
the host galaxies (e.g., Bloom et al. 2001). Observed GRB 
hosts are small-mass galaxies, often thought to be going 
through vigorous star formation phase. 

There are no reliable GRB host mass estimates, and thus 
we have calculated models for a range of galaxy masses. 
Our standard model calculations were repeated for a num¬ 
ber of different evolutionary models to assess the robust¬ 
ness of our results. We have found that the NS-BH merg¬ 
ers take place mainly outside of their host galaxies, and 
thus are inconsistent with the observed locations of GRBs 
around hosts. Some WD-BH binaries may merge outside 
the star formation regions of their host galaxies. However, 
the distribution of the WD-BH merger sites around their 
host galaxies is consistent with the observed distribution 
of GRB offsets from the centers of galaxies. Thus one can 
not reject the WD-BH mergers purely on the basis of com¬ 
parison with the observed offsets. However, if one addi¬ 
tionally requires that the mergers should take place in the 
proximity of star forming regions in galaxies then the WD- 
BH mergers can be rejected as potential GRB candidates. 
Merger sites of WD-NS, He-NS, He-BH, and NS-NS trace 
the star formation regions of the hosts, for all the cases of 
the host mass and size considered here, and independently 
of the adopted population synthesis model. We conclude 
that these types of binaries may be responsible at least for 
a part of observed GRBs. 

GRBs form a very nonuniform group of events, with 
different outburst times, very different light curves and 
observed energies. Thus, there is a possibility that GRBs 
originate in more than one type of progenitor. Locations of 
GRBs with respect to host galaxies has so far been mea¬ 
sured only for long GRBs. There is a growing evidence 
that these GRBs are related to collapsing massive stars: 
collapsars. However, our results show that several types 
of binary system progenitors cannot be rejected purely on 
the basis of their merger site distribution. Additionally, if 
binaries were responsible for only a part of the observed 
GRBs, we also can not exclude them purely on the basis 
of their expected coalescence rates. 

Because of the expected short duration times, NS-NS 
and NS-BH mergers are the primary candidates as short 
burst progenitors. These two populations exhibit very dif¬ 
ferent distributions of merger sites. Mergers of NS-NS sys¬ 
tems take place predominantly within hosts, to the con¬ 
trary of what was so far believed, provided that CE phases 
initiated by low-mass helium stars do not always lead to 
binary component mergers, assumption which yet have to 
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be tested by detailed hydrodynamical calculations. On the 
other hand, a significant fraction of NS-BH systems merge 
outside of their host galaxies. At some point in the future 
afterglows from short GRBs will be observed and their 
locations with respect to host galaxies will be measured, 
and then such calculations may provide a useful tool to 
distinguish between these two progenitor models (Perna 
& Belczynski 2001). Future and current space missions 
like HETE-II, INTEGRAL, GLAST or SWIFT will hope¬ 
fully measure precise positions of a large number of bursts 
even of the short duration and settle down the issue of 
GRB progenitors. 
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Fig. 1. White dwarf mass distributions in WD-NS (solid line) and WD-BH systems (dashed line) for our standard 
evolutionary model. Distributions are normalized to the total number of binary GRB progenitors (100800) formed out of 
IVtot = 3 x 10 7 primordial binaries. We require that a WD mass exceeds 0.9 M 0 for a given system to be classified as a 
potential GRB progenitor. We also study models in which the minimum mass is 0.7 and 1.1 M 0 . 


Table 1 

Characteristic Timescales of GRBs Candidates (Myr) 


Type 

^evol 

A^evol 

t b 

t'merg 

^^merg 

WD-NS 

26.6 

20.0-35.7 

6.8 

0.014-1238 

WD-BH 

25.6 

20.1-34.6 

96.9 

16.1-1831 

He-NS 

10.0 

7.04-11.1 



He-BH 

7.9 

5.62-9.96 



NS-NS 

18.5 

10.7-27.8 

0.7 

0.017-390 

NS-BH 

7.7 

5.92-17.4 

534.6 

1.68-5170 


distribution median of evolutionary time delay 
b distribution median of merger time delay 
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Table 2 

Population Synthesis Model Assumptions 


Model Description 

A standard model described in § 2 

Bl-13 zero kicks, single Maxwellian with 

cr = 10,20,30,40, 50,100, 200, 300,400, 500, 600km s" 1 , 

“Paczynski” kicks with a = 600kms _1 
C no hyper-critical accretion onto NS/BH in CEs 

Dl-2 maximum NS mass: M maXi ns = 2,1.5M 0 
El—3 ccce x A = 0.1,0.5,2 

FI -2 mass fraction accreted: f a = 0.1,1 
Gl-2 wind changed by / w i n d = 0.5, 2 
H Convective Helium giants: M conv = 4.0 M 0 

I burst-like star formation history 

J primary mass: oc M~ 2 - 35 

Kl-2 binary fraction: /b; = 0.25,075 

LI -2 angular momentum of material lost in MT: j = 0.5, 2.0 
Ml-2 initial mass ratio distribution: <f>(g) oc q~ 2 7 ,q 3 
N no helium giant radial evolution 

O partial fall back for 5.0 < Mqo < 14.0 M 0 

Pl-2 minimum Helium core mass in He-NS/BH mergers A7 m j n H e = 4, 8 M 0 
Rl-2 minimum WD mass in WD-NS/BH mergers: M m j n ,WD = 0.7,1.1 M 0 


10- 1 


IQ-2 


^ 10- 3 


10- 4 


10-5 


Fig. 2.— Helium core mass distributions in He-NS (solid line) and He-BH mergers (dashed line) for our standard 
evolutionary model. Distributions are normalized to the total number of binary CRB progenitors (100800) formed out of 
A^tot = 3 x 10' primordial binaries. We require that a He core mass exceeds 6M 0 for a given system to be classified as 
a potential CRB progenitor. We also study models in which the minimum mass is 4 and 8 Mq. 
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Table 3 

Galactic Binary GRB Progenitors Coalescence Rates (Myr -1 ) 


Model 3, 

WD-NS 

WD-BH 

He-NS 

He-BH 

NS-NS 

NS-BH 

A 

4.6 

2.4 

9.7 

23.5 

52.7 

8.1 

B1 

46.3 

13.4 

20.9 

64.2 

292.4 

18.2 

B2 

50.9 

12.9 

21.4 

62.9 

299.6 

19.4 

B3 

48.7 

13.6 

20.8 

63.7 

302.2 

19.6 

B4 

44.6 

12.2 

20.7 

66.5 

285.2 

19.1 

B5 

38.2 

11.3 

22.8 

67.2 

251.0 

19.5 

B6 

32.2 

10.3 

19.9 

64.0 

226.8 

16.4 

B7 

13.4 

5.3 

15.2 

48.9 

128.1 

14.6 

B8 

4.8 

2.6 

9.9 

23.5 

57.5 

10.1 

B9 

1.9 

0.9 

8.9 

12.8 

33.2 

5.7 

BIO 

0.8 

0.9 

6.9 

9.7 

18.2 

3.7 

Bll 

0.4 

0.4 

6.2 

7.6 

12.5 

2.1 

B12 

0.4 

0.4 

4.6 

5.8 

8.2 

1.5 

B13 

12.2 

4.0 

12.3 

29.8 

91.0 

10.3 

C 

0.4 

1.7 

33.3 

12.7 

43.2 

5.6 

D1 

104.8 

7.7 

1.8 

33.8 

33.6 

23.3 

D2 

114.7 

22.1 

0.1 

32.4 

9.1 

36.2 

El 

0.03 

0.2 

0.5 

91.6 

2.5 

4.7 

E2 

1.7 

0.3 

8.5 

47.8 

23.5 

6.3 

E3 

5.4 

6.0 

4.6 

8.1 

109.0 

8.7 

FI 

4.2 

2.1 

2.3 

14.5 

22.1 

9.3 

F2 

6.5 

11.1 

8.2 

4.5 

54.3 

8.6 

G1 

5.7 

5.8 

7.2 

20.3 

43.9 

14.2 

G2 

4.8 

0.6 

19.7 

15.1 

94.8 

1.3 

H 

4.7 

2.0 

8.2 

24.3 

37.9 

7.8 

I 

4.3 

3.5 

9.7 

23.9 

54.5 

10.0 

J 

4.8 

3.8 

12.6 

34.8 

58.1 

12.8 

K1 

1.9 

1.0 

4.1 

9.9 

22.5 

3.4 

K2 

7.8 

4.0 

16.3 

39.6 

90.2 

13.5 

LI 

6.0 

3.6 

6.9 

8.4 

78.9 

9.2 

L2 

4.3 

2.2 

6.6 

21.3 

12.0 

6.2 

Ml 

0.9 

4.3 

1.2 

5.9 

6.2 

4.0 

M2 

5.8 

0.2 

17.4 

22.9 

114.2 

8.4 

N 

7.5 

4.0 

8.7 

22.3 

34.4 

10.7 

0 

4.2 

1.6 

10.0 

25.0 

51.9 

5.7 

PI 

4.6 

2.4 

73.6 

33.3 

52.7 

8.1 

P2 

4.6 

2.4 

0.9 

10.7 

52.7 

8.1 

R1 

5.6 

3.5 

9.7 

23.5 

52.7 

8.1 

R2 

2.9 

1.4 

9.7 

23.5 

52.7 

8.1 


for definition of models see Table || 
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Table 4 

Location of GRB afterglows in relation to their host galaxies a 


GRB 

redsliift 

Offset AO 

Rprojected [kpc] 

Comments 

970228 

0.695 

0.426 ± 0.034” 

3.266 ±0.259 

edge of host 

970508 

0.835 

0.011 ±0.011” 

0.091 ±0.090 

host center 

970828 

0.958 

0.474 ±0.507” 

4.047 ±4.326 

edge/outside 

971214 

3.418 

0.139 ±0.070” 

1.105 ±0.557 

inside host 

980326 

~ 1 

0.130 ±0.068” 


edge/outside? 

980329 

£3.5 

0.037 ± 0.049” 


inside host 

980425 

0.008 

12.550 ± 0.052” 

2.337 ±0.010 

inside host 

980519 


1.101 ±0.100” 


inside host 

980613 

1.096 

0.089 ± 0.076” 

0.782 ±0.666 

??? 

980703 

0.966 

0.040 ±0.015” 

0.038 ±0.128 

inside host b 

981226 


0.749 ±0.328” 


??? 

990123 

1.600 

0.669 ± 0.003” 

6.105 ±0.027 

edge of host 

990308 


1.042 ± 0.357” 


??? 

990506 

1.310 

0.297 ± 0.459” 

2.680 ±4.144 

??? 

990510 

1.619 

0.066 ± 0.009” 

0.600 ±0.084 

edge of host 

990705 

0.840 

0.872 ± 0.046” 

7.165 ±0.380 

inside host 

990712 

0.434 

0.049 ±0.080” 

0.301 ±0.486 

inside host 

991208 

0.706 

0.196 ±0.097” 

1.513 ±0.750 

edge? 

991216 

1.020 

0.359 ± 0.032” 

3.107 ±0.280 

inside? 

000301C 

2.030 

0.069 ±0.007” 

0.622 ±0.063 

inside? 

000418 

1.118 

0.023 ± 0.064” 

0.202 ±0.564 

host center 

000926 

2.066 

1.5 ± 0.5” 

13.43 ±4.5 

edge/inside c 

010222 

1.477 

0.05 ± 0.05” 

0.45 ± 0.45 

inside host d 


a all data from Bloom et al. (2001), but: 

b Berger et al. (2001) 

c Fynbo et al. (2001) 

d Jha et al. 2001; Fruchter et al. 2001 
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Table 5 

KS TEST COMPARISON RESULTS BETWEEN MODELS AND OBSERVED OFFSETS 3 


Galaxy mass: 

Mmw 


0.1 Mmw 

O.OIMmw 

O.OOIMmw 

WD-NS mergers 

Standard model 

2.42 x 10" 4 


3.50 x 10“ 2 

5.87 x 10~ 4 

9.48 x 10” 1 

maximal model b 

2.94 x 1CT 4 

(B7) 

4.41 x 10~ 2 (B8) 

5.29 x 10" 1 (LI) 

3.26 x 10" 1 (F2) 

minimal model b 

1.20 x 10” 4 

(F2) 

2.53 x 10” 2 (F2) 

8.04 x 10" 2 (B12) 

6.54 x 10" 4 (B12) 

WD-BH mergers 






Standard model 

2.24 x 10" 4 


3.49 x 10" 2 

5.87 x 10” 1 

9.48 x 10" 1 

maximal model 

2.86 x 10~ 4 

(B6) 

3.92 x 10" 2 (B2) 

6.85 x 10 _1 (E2) 

9.57 x 10 _1 (J) 

minimal model 

1.01 x 10“ 4 

(E3) 

1.80 x 10~ 2 (E3) 

8.23 x 10" 2 (G2) 

6.65 x 10 -4 (G2) 

He-NS mergers 






Standard model 

2.70 x 10" 4 


3.61 x 10" 2 

8.27 x 10~ 2 

6.71 x 10~ 4 

maximal model 

3.57 x 10~ 4 

(B7) 

4.39 x 10 -2 (C) 

8.32 x 10 -2 (J) 

7.12 x 10" 4 (FI) 

minimal model 

2.22 x 10~ 4 

(E3) 

3.57 x lO’ 2 (FI) 

6.95 x 10 -2 (B6) 

6.21 x 10" 4 (Bl) 

He-BH mergers 






Standard model 

2.84 x 10" 4 


3.81 x 10~ 2 

8.15 x 10 -2 

7.07 x 10~ 4 

maximal model 

3.21 x 10“ 4 

(FI) 

3.99 x 10 -2 (G2) 

1.00 x 10 _1 (I) 

7.98 x 10 -4 (L2) 

minimal model 

1.71 x 10” 4 

(Bl) 

3.15 x 10” 2 (B4) 

7.87 x 10 -2 (B7) 

6.28 x 10" 4 (Bl) 

NS-NS mergers 






Standard model 

1.84 x 10" 4 


2.90 x 10" 2 

3.39 x 10" 1 

1.04 x 10~ 2 

maximal model 

2.44 x 10~ 4 

(B2) 

3.63 x 10 -2 (Bl) 

5.72 x 10" 1 (M2) 

2.86 x 10" 1 (F2) 

minimal model 

1.20 x 10 -4 

(F2) 

2.14 x 10~ 2 (F2) 

1.62 x 10 _1 (Bl) 

2.61 x 10" 3 (Bl) 

NS-BH mergers 






Standard model 

3.13 x 10~ b 


4.69 x 10" 4 

1.25 x 10~ b 

8.44 x 10" s 

maximal model 

2.72 x 10~ 4 

(D2) 

3.16 x 10" 2 (D2) 

4.07 x 10" 1 (Dl) 

2.19 x 10" 2 (Dl) 

minimal model 

8.37 x 10“ 6 

(0) 

1.80 x 10“ 6 (0) 

9.78 x 10 _1 ° (0) 

1.09 x 10 _1 ° (0) 


a We list the probabilities that the observed offsets distribution has been drawn from the theoretical one 
b corresponding models are given in parenthesis 
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Fig. 3.— Distributions of evolutionary (broken lines) and merger (solid lines) times for GRB binary candidates calculated 
in our standard evolutionary scenario. Distributions are normalized to the total number of binary GRB progenitors 
(100800) formed out of TVtot = 3 x 10 7 primordial binaries. Note that every panel has diffrent vertical and horizontal 
scales. 
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Fig. 4.— The dependence of galactic GRB progenitor coalescence rates on the assumed natal kick velocity distribution. 
Lines connect rates for models B1-B12 and the horizontal scale shows the width of Maxwellian kick distribution of a given 
model. Triangles mark rates of our standard model (A) and squares he model with ‘Paczynski” kick distribution (B13). 
The width of kick velocity distribution scale is irrelevant for these two models, and they were placed in horizontal axis to 
approximately match the rates obtained with single Maxwellian kick velocity distribution. 



Fig. 5.— Star formation history rates used in this work. The thin line is based on lower limits from Madau et al. (1996), 
while the thick line represents approximately the rate of Rowan-Robinson (1999). 






20 




Fig. 6.— Cumulative event rates of different GRB progenitor types as a function of redshift for our standard evolutionary 
model. From top to bottom curves correspond to: NS-NS (dotted-short dashed line), He-BH (long dashed line), He-NS 
(short dashed line), NS-BH (solid line), WD-NS (doted-long dashed line) and WD-BH mergers (dotted line). The left 
panel shows the case with assumed star formation rate history of Rowan-Robinson (1999), while the right panel to that 
of Madau et al. (1996). For all calculations flat cosmology model was used, with fl m = 0.3 and = 0.7. 








rate [yr *] rate [yr 1 \ rate [yr 



Fig. 7. Cumulative event rates of different GRB progenitor types as a function of redshift for all our different models. 
All rates were calculated with assumed star formation rate history of Rowan-Robinson (1999). For all calculations flat 
cosmology model was used, with f2 m = 0.3 and Oa = 0.7. 
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Fig. 8.— Cumulative distributions of neutron star black hole binaries merger sites around different mass galaxies (solid 
line) for our standard evolutionary scenario (model A). The initial distribution of primordial binary population within 
the galaxy is shown with the dashed line. 
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Fig. 9.— Cumulative distributions of white dwarf black hole binaries merger sites around different mass galaxies (solid 
line) for our standard evolutionary scenario (model A). The initial distribution of primordial binary population within 
the galaxy is shown with the dashed line. 
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Fig. 10. - Cumulative distributions of several types of GRB progenitors around galaxies for our standard evolutionary 
scenario (model A). The case of white dwarf neutron star mergers is illustrated in the top panel with to extreme cases: a 
Milky Way like galaxy (Left), and small galaxy with the mass O.OOIMmw- We present the distributions of Helium star 
mergers in the middle panel for the case of a small galaxy only. The lower panel contains the plots with the distributions 
of double neutron star mergers around a Milky Way like galaxy (Left), and around a small dwarf galaxy with the mass 
O.OOIMmw- The initial distribution of primordial binary population within the given galaxy is shown with the dashed 
line. 
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Fig. 11. - Cumulative distributions of WD-NS and WD-BH merger sites for two extreme galaxy masses and for different 
evolutionary models. All models are shown with the dotted lines, except the most extreme ones: model F2 - short dashed 
line, N - doted long dashed line, E3 - short long dashed line, LI long dashed line. The initial distribution of primordial 
binary population within the galaxy is shown with the solid line. 
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Fig. 12.— Cumulative distributions of NS-NS and NS-BH merger sites for two extreme galaxy masses and for different 
evolutionary models. All models are shown with the dotted lines, except the most extreme ones: model F2 - short dashed 
line, N doted long dashed line, E3 - short long dashed line, D1 dotted short dashed line, D2 long dashed line. The 
initial distribution of primordial binary population within the galaxy is shown with the solid line. 








